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ABSTRACT 


A previously developed program, which includes all 
electronic interactions thought to be important, does not 
correctly predict the value of electron mobility in mercury 
cadmium telluride particularly near room temperature. Part of 
the reason for this discrepancy is thought to be the way 
screening is handled. It seems likely that there are a number 
of contributors to errors in the calculation. The objective 
of this work is to survey the calculation, locate reasons for 
differences between experiment and calculation, and suggest 
improvements . 
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1 . 


Introduction 


Many relevant properties of and references about 
Hg, Cd Te have already been given (Patterson (1987), see also 
Appendix 1) . 

The calculation of the transport properties of narrow 
band gap semiconductors (of which mercury cadmium telluride is 
an important example) have met with considerable success 
(Lehoczky et al. (1982), Szymanska et al. (1978), and 
Chattopadhyay and Nag (1974)). Nevertheless it is not unusual 
to find calculations (near room temperature) differing from 
experiment by a factor of order two. 

The original idea of this project was that inadequate 
treatment of screening was the cause of discrepancies between 
calculation and experiment. However, in a calculation as 
complex as the mobility of electrons in a narrow band gap 
semiconductor (see Appendices 2 and 3) even the accuracy of 
the present calculations whould have to be regarded as good 
(Lehoczky et al. (1980)). 

An easy way to test the importance of screening is to run 
the present codes with different values of the screening 
constants. As we will see such calculations lead one to feel 
there is a limit to the accuracy with which the mobility can 
be evaluated, and certainly a limit to how much the accuracy 
can be improved by improving the screening. 
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2 . 


Objectives 


a. Determine the effect on the mobility of electrons in 
Hg, Cd Te of variations in relevant parameters. Pay particu- 
lar Ittintion to screening parameters. 

b. Survey the theory on which the calculation is based 
to try to pinpoint areas that can be improved. 
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3. 


Calculations 


We start by observing that in the approximation con- 
sidered, the effects of screening are determined by the 
screening lengths. There exist reasons for considering other 
values of these lengths beside the conventional one 

(Patterson, (1987), appendices 3 and 5). Thus we are 
motivated to consider the results of calculations in which the 
screening length is varied for fixed values of the other 
parameters . 

We also show how variations in screening affect the 
results for different values of the Cd concentration (x) . All 
our results are based on a program of Dr. S.L. Lehoczky as 
rewritten by Ernestine Cothran. The program considers 
electron interaction effects due to longitudinal optic 
phonons, ionized impurities, holes, acoustic phonons, composi- 
tional disorder and neutral defects. 

In what follows, we present our calculations in terms of 
hypotheses and calculational tests of the hypotheses. 
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Hypothesis 


Although screening is important sizable changes in 
screening are not sufficient to get agreement with experiment. 


TEST (iA 
vector's 


is the mobility, k T „ represents all Fermi-Thomas wave 
and T is the temperature) . 



( , £> 7 -f=- J )( — . / i 3 ) — 
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at T 


at T 


14 -3 

300K, ^ 5.6 x 10 cm (the donor concentration) 

f E = .742, 

100K, n p = 5.6 x 10 14 cm -3 


f £ ~ .665. 

This is further shown by Figs, la and lb where one sees that 
it is necessary (for the example shown) to use .05 of the 
nominal values for j\.~ in order to obtain near agreement with 
experiment. t 
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Hypothesis 


as a 


Test 
at x 


at x : 


The mobility varies monotonically (very roughly linearly) 
function of inverse screening length (k TF ) . 
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Hypothesis 


The more carriers, the more screening we should have. 
Thus, the more the calculation should depend on variations in 
the screening length. Three ways the number of carriers can 
be varied are by varying the Hg concentration, varying the 
temperature T or varying the donor concentration. 

Tests (These imply the above hypothesis is naive.) 

at T= 300K, 5.6 x 10 14 cm" 3 


/, l ( „ )(=J33j 


at x= .183, 17 = 5.6 x 10 14 cm 3 
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At a fixed temperature, changes due to screening are more 
important for smaller band gaps. This is shown clearly in the 
Figs. 2a and 2b. 

A simple expression for the energy gap as a function of 
temperature and x is given by Dornhaus and Nimtz (1983) who 
also list more accurate expressions. The expression is only 
supposed to be used for .135 < x < .203. In the table below 
the energy gap is evaluated from a slightly better expression 
in the program. 

E g = -.3 + 5 x 10" 4 T +(1.91 - 10 -3 T)*x in eV. 

We find at T = 300K, 

E = -.15 + 1.61X. 

g 

If x is the value of x for which the energy gap 
vanishes, ^e find for T = 100 that x q = .138. 
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FIG 2. (A) ENERGY GAP E g (eV) AND CHANGE IN MOBILITY aM ( Cm 2 / Vs ) '* 

(B) FRACTIONAL CHANGE IN MOBILITY VS ELECTRON CONCENTRATION. 


4 . Conclusions and Recommendations 


a. Screening is important but it is not the whole 
answer. There are undoubtedly several small corrections. 

b. In a calculation as complex as this one, see Figs. 3a 
and 3b and the appendices, the existing calculation (Lehoczky 
(1980), (1982)) is already good. See also Nag (1980) for 
comparison. 

c. Small corrections and rewriting the code for a more 
modern computer are realistic goals for improvement. The 
program is presently implemented in basic for an HP 9845B and 
it runs rather slowly. 
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Appendix 1 

Brief Properties of Ha 1 _ ^Cd^Te: 

1. Pseudo Binary II-VI Compound with structure isomorphic to 
zincblende. 

2. The bandgap is tunable form 0 to 1.605 eV as x goes from 
about .17 (at low temperature) to 1.0. The bandgap also 
depends on temperature. Hg Cd Te forms a continuous range 
of solid solutions between the seMmetal HgTe and the semicon- 
ductor CdTe. 

3. Can use as an infrared detector at liquid nitrogen 
temperature in the wavelength range 8-12^ - an atmospheric 
window. A higher operating temperature than alternative 
materials, high detectivity, fast response, high sensitivity, 
IC compatible and low power. 

4. Candidate for microgravity crystal growth. 

5. Band structure involves mixing of unperturbed valence and 
conduction band wave functions as derived by the Kane theory. 
Non parabolic bands. 

6. Typical carriers have small effective mass (about 10 -2 * 
free electron mass) which implies large mobility. The 
effective mass also has a strong dependence on concentration 
and temperature. 

7 . Polar scattering by LO lattice vibration modes dominate 
above 77K (or so) and the relaxation time approximation is not 
valid. The scattering is inelastic since the phonon energy is 
comparable to the electron energy. 

8. Narrow bandgap semiconductors imply thermally excited 
electrons and hence screening. There may also be 
antiscreening effects when the phonons have higher phase 
velocity than the Fermi velocity of electrons. 

Dornhaus and Nimtz (1983) list many more properties. 
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Appendix 2 

Boltzmann Transport Equation and Solution 


We just have a few comments to make here as the Boltzmann 
Transport equation is well reviewed in several texts (Blatt 
(1968), Wilson (1958)). 


If f is the distribution function (with f its value in 
equilibrium), E the energy, and if terms of ° order flF 2, are 
neglected where 


S- * £ 


( k ) 


i£o 

oe 


(ftt. i) 


then with electric fields £ 
the x-direction 


^x 



■r 


5-V 


and temperature (T) 



gradients in 




a/T', 


where E p is the Fermi energy v is the velocity, k the wave 
vector, /<?/ the magnitude of electronic charge, and 
V(^,F' ) is determined by the electron scattering mechanism 
(it_^ is proportional to the probability of scattering from 
to it * ) . 


Since some of the scattering mechanisms are inelastic a 
variational principle is used to solve for &(hj . 

The work is detailed, but well described ((Wilson (1958) 
Lehoczky (1980)). From f, the current density J 
(with 2l=o ) is X 



( 42 . 


and the mobility is M - Jx/f h/et?*.), 
where y) is the electron concentration. 
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Appendix 3 

Electron Scattering Mechanisms 


The details here are too long to write down. Fortunately 
they have already been well described (Lehoczky (1980) , 
Lehoczky (1974), Nelson (1978), Dubowski (1981), Ehrenreich 
(1959), Zawadzki (1974), and Szymanska (1978)). Special 
mention should be made of Ehrenreich who wrote many of the key 
papers in this area. They are referenced in the above 
references. 

a. Electon-phonon Interaction in Polar Solids 


The only strong interactions are with longitudinal 
optic modes and we will limit discussion to long wave length 
modes near the center of the Brillouin zone with frequency^ . 
Several approximations then follow for which the literature 
may be consulted. 




M 

f ( ) 



( tea, 0 


The a's are the phonon annihilation and creation operates and 
the c's are the electron creation and annihilation operators. 
M is an apropriately chosen constant. 


The screening is included by using 

function 



the 


dielectric 


(fl 3 6>. <?J 


and a "two mode" theory is necessary (with frequencies close 
to those of CdTe and HgTe) . £ is the high frequency 

dielectric constant. v< 
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Further Comment Screening and Antiscreenina 


Our main point here is to point out that antiscreening at 
small wave vectors can account for increased interaction and 
hence reduced scattering. This can occur when phonons with 
small wave vector also have high phase velocity which can 
occur for optical modes. As Ehrenreich (1959) has shown, the 
screened potential can actually be greater than the unscreened 
one. 


b. Charged Impurity and the Friedel Sum Rule 


This is also referred to as ionized impurity 
scattering ( ii ) : 


H ,7 = Hi Sp c *r C * v > 


Hi; 

1 Hi 


c r 

where 

F r, r 

with 

Vi? (r) 


XL 


(r) Vf , 


( flU'l) 


( /lib. ?) 


in Bloch's theorem and is often not included. V(q) is the 
Fourier transform of the Coulomb potential. Screening is 
included by £(q) and 


- 


* U 


f + (C) z 


'FT 


Ft 



where is the low frequency dielectric constant. 


Screening was more fully discussed in the first report 
(Patterson, 1987) . Briefly, if a charge is placed in an 
electron gas its potential is better described by cf'* r /r' 
than */r because the electrons in the gas surround the change 
and "screen" it. The screening length is defined as 

. The process is best described by a wave vector depen- 
dent dielectric constant, as above. 
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Further Comments: Multiple Scattering and Other Corrections 
to Calculations 

Chattopadhyay and Queisser (1981) have discussed cor- 
rections due to multiple scattering and impurity dressing. 
The usual assumption with regard to ionized impurity scatter- 
ing is that each impurity scatters independently. Multiple 
scattering adds corrections due to the electron interacting 
with more than one ion simultaneously. There is also a 
correction due to the fact that the impurities perturb the 
electron energies and wave functions. This is called the 
dressing correction and it can amount to 50% or so of the 
multiple scattering correction. However the corrections to 
ionized impurity scattering are not very important compared to 
other scattering mechanisms in Hg Cd Te at temperatures of 
interest here. The dressing ana multiple scattering cor- 
rections to ionized impurity scattering are more important at 
higher temperatures and lower impurity concentrations. 

Corrections to the Born approximation (in a calculation 
of ionized impurity scattering) can also be important at lower 
temperature but second order Born corrections are probably 
only useful if they are small. 

An expression for the screening length L = PC 1 is given in 
appendix 5. It is possible to derive this directly or by use 
of the first Born approximation and the Friedel sum rule 
involving phase shifts fa (see e.g. Joachain, 1975) . The 
Friedel sum rule holds generally and is supposed to be exact. 
It holds even for interacting electrons in a periodic sysem. 
The essential physical contents of the Friedel sum rule is the 
requirement of overall electrical neutrality about the charged 
ion. For electrons in a semiconductor, the Friedel sum rule 
takes the form, in the usual rotation; 

2 = (z/jt) 2Jj, f 2§rfe> ■ V 


c. Electron-Hole 


When the heavy hole mass is large compared to the 
electron mass the hole can be considered to be stationary 
during scattering. This would make electron-hole scattering 
theory formally identical to ionized impurity scattering. 

However we still must include non parabolity and 
admixture of p wave functions into the conduction band. For 
details see the references. 
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This is also a screened interaction, 
dielectric function has been used 

yZ / , | 2 / > 4* 1 * 

) r ( L r ) 


The following 




0< 


e 

Fr 


-t 


(&) 


( A3C.t) 


d. Electron-Acoustic phonon 

For semiconductors this is often treated by using the 
idea of a deformation potential (Seeger 1985, Madelung, 
1978) . If E c is the bottom of the conduction band then if c. 
is the energy of the electron, we can obviously write 


£ = (e- e. ) i-F 




where (F-^o ) is treated as the kinetic energy and is 

viewed as the potential energy. One can think of varying 

with long wavelength displacements ( ) and this variation 

produces the deformation potential and thus the electron 
lattice coupling. At long wavelengths, Umklapp processes can 
be neglected. We write 


U. 


/ - pk 


df 0 

57 


AV - V V- <P ^ F* {A 34*2) 
7> y o > 


for a single wave with wave vector <5 . Further generaliza- 

tion leads to putting this in the usual form (see references) . 
In the above V is the volume. 

In the actual calculations, both longitudinal (L) and 
transverse (T) scattering is included by the use of three 
deformation potentials (1L, 2T) . 

e. Compositional Disorder and Neutral Defect 

Compositional Disorder 

This is calculated by considering a random distribu- 
tion of scattering centers consisting of square wells. 



XXII-20 




The square wells are chosen to have a size equivalent to a 
unit cell and a depth (E.) determined by the differences in 
band gaps between CdTe ancr HgTe . 

If Na is the number of unit cells per unit volume 

£/ ) xo-x) ej- 

'Jfs. k.,TT > (fl3e. /) 

c? 

(see Lehoczky et al., (1980)). 

The dependence on x(l-x) is known as Nordheim's rule. For a 
given x, x(l-x) measures the disorder assumed to be present 
(maximum) . 

Neutral Defect Scattering 

This is also described as scattering by a random distri- 
bution of square-well potentials. It has been neglected in 
this calculation, but one assumes 

— - ) c*Xfl-*)NE V . (Alt. 2) 

(see Lehoczky et al., (1980)), 

where N is the density of neutral defects, E is their strength 
and V is the volume of the cell. 
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Appendix 4 

The Polaron Constant 


The Polaron 
notation, as 


constant °< is defined, 







in the usual 


(fHj) 


It is a measure of the strength of the electron-lattice 
interaction and is dimensionless. For the phonon field and 
electrons to interact weakly we need x <1. For the Frohlich 


interaction 

term to 

be valid we 

also 

need x <1. 

To get 

an estimate of °C we 

use 

the following values: 

T(cO) 

= 15 

9 

lei 

= 4.8 x 10“ 10 esu 

ft o) 

= 20 

t 

fr 

= 1.05 x 10~ 27 erg-s 


— . 03m 

9 

hi 

= 9.1 x 10 -28 gm 


= 2.4 X 

-14 

10 ergs 




Thus 

= .092. 




I 
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Appendix 5 

Band Structure and Effective Masses 


We briefly summarize Kane (1957) theory here. We give 
some idea of the band structure near the top of the valence 
band and bottom of the conduction band. A simple schematic 
diagram of the band structure has been given by Seeger (1985), 
see also Lehoczky (1974) . The s wave functions of the conduc- 
tion band are found to mix with the p wave functions of the 
valence band by perturbation theory (for 7z £0 ) . 

By using Bloch's theorem the wave function is given by 

e ‘" r - ffs-.O 


in the usual notation. The Schrodinger wave equation can be 
written: 

('u.i-H') u„ r (?) = e,' u kF (?) t C#s.z) 


where 


- 

A/’ = 


r 


t | l(r) 

i. t- p + (Vt'xp)-r', C^ -yJ 

ton 1 ° ' 


and a \7 V k k 


term has been neglected. 

£ = £ — * /(2ho) . 


O'. 


(/)£.£) 


The Pauli operation is <7" . 


To solve for the required energies, the eigenvalues of an 
8x8 matrix must be determined. 
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ORIGINAL PAGE IS 
OF POOR QUALITY 


In Kane's notation we define: 

p = -Ul)< s //*/*> , 


(4 ST. 6 ) 


and 


Sir,' 


A ^ ~ P A/^>, (**?) 




DP 

'y 


where the matrix elements involve s and p ( / / Ty ^ / ?> ) 

wave functions. If A?>fcP, as seems to be appropriate for 

narrow band gap semiconductors, we get for the energy as 
measured from the bottom of the conduction band (neglecting a 
small term ti* k l /a m ): 

£■ = - +- \[(£j/ 3 f-r , (VAp 


where is the energy gap between the conduction and the 

valence band. 


Depending on the usage, various effective masses can be 
defined for spherical energy surfaces as above. We assume 
spherical energy surfaces and consistant with above 
(Chattopadhyay and Nag, (1974)): 

E ( 

where m 0 is the mass at the bottom of the conduction band 



9 A/ P*- 



Energy Effective Mass 

pf, 



i? E 


CfiS’./o) 


Energy Effective Mass 2, 

m £z = 


t* / — r‘ 


(rfs: j/j 
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Momentum Effective Mass, 


m r =■ Pk ( ' ' fjs-./z) 

( <j- = £ k and V~ and ^ are parallel) . 

It is relatively easy to show 

?»>,/*,„ -/ +■ , (/fz.aj 

'"/>/*» - / -t se/, ^s > (/?$./•/) 

2 . _/ 

Y"£?/m 0 - (l r 2£ /e ) f (Jsr./sj 

all three of these agree at the bottom of the conduction band 
( £T ~ o } J? — o ) • 

Which one of these is appropriate depends on the usage. 


The conventional effective mass as defined, for example, 
by Kittel (1986, p 194) is togz. • The effective mass equation 
may be used as a starting point for discussing scattering in 
the Born approximation in solids. It would appear from a 
development of that equation (Madelung (1978) , p. 79) that 
energy effective mass 1 could be the appropriate effective 
mass to use. For HgCdTe, Wp is typically about 10 of the 
mass of free electrons and has a strong dependence on carrier 
concentration and temperature. 

By a standard development in the first Born approxima- 

Patterson 




parabolic 

momentum 

properly 


r/J 


tion, the screening length comes out (see e.g. 
(1987)), 


L 




7T% 




However, for this result to be applicable to non 
energy bands, the effective mass should be the 
effective mass so that the density of states (P(£) ) 
enters via 


p(£)Je =■ 


7T 


)e 


• d_ k 


ds s 


J£ 


TT^-k 2 - 
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Since ft) 


~k (top rAi e ) 


, (see above) 


£7 


we find: 


L 


- z 
to e < 





Then since 


top > to 0 , 




we have 


z. 


-z. 

^/T/ 



(/?£■)*)) 


Hence a crude approach using the Born approximation could 
lead to a larger screening length, hence more scattering, 
hence less mobility. 

Further Comment: 


We also need to point 
an infrared (IR) detector. 


out the usefulness of Hg Cd as 
Since 1-x x 


F = ~kco ~ he /* j 


it is easy to show that 


e„ (et‘) = / , 


2u) 


An important atmospheric window for IR is in the wavelength 
range 8-12^ . At T=77K the energy gap is approximately 
(Dornhaus and Nimtz, 1983) 

£■; - ~ .2<* r/.f?X . (AS\2t) 

a 

We find that .16 £ x £ . 19 will put us in the 8-12/ c window. 
There are other possibilities. ^ 
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Appendix 6 
Miscellany 


a. Wiener Lattice 

In an electron gas in which the kinetic energy is much 
less important than the potential energy, a localization or 
"crystallization" of the electrons can occur. It appears that 
this phenomena actually occurs at x = .24 in Hg Cd Te at 

millikelvin temperature in a magnetic field (Rosenbaum <lt al. 
(1985)). 

b. Layered Structures 

One can make a superlattice consisting of alternating 
layers of CdTe and HgTe by use of molecular beam epitaxy * 
Such materials can be used as infrared detectors with proper- 
ties that may be superior to Hg Cd Te alloys (Berroir and 
Voos (1987)). 

c. Hg 1 _ x Zn x Te 

As this report reflects, Hg Cd Te is the most important 
intrinsic semiconductor infratec^ detector material. 
Hg. Zn Te has been considered as a possible competitor, 
because* it may be more structurally stable. The use of 
Hg, Zn Te is still quite new and much effort is being devoted 
to 1-3 Sri^tal growth techniques (Jozwikowski and Rogalski 
(1988) ) . 
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